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Abstract 

We review the different frameworks in which Galileon scalar fields have been seen 
to emerge such an in DGP, New Massive Gravity and Ghost-free massive Gravity 
and emphasize their relation with the Lovelock invariant in braneworld models. The 
existence of a non-renormalization theorem for Galileon scalar fields makes them 
especially attractive candidates for inflation as well as for late-time acceleration. In 
particular we review the self-accelerating and degravitating branches of solutions 
present in Galileon models when arising from Massive Gravity and discuss their 
phenomenology. 

1 Introduction 

The late-time acceleration of the Universe remains undeniably one of the most in- 
triguing aspects of modern Cosmology, at the interface between General Relativity 
and particle physics. Over the past decade a multitude of new models have been 
developed all across the community to tackle this fascinating issue, see Refs. [1, 2, 3] 
for different reviews as well as Refs. [4, 5, 6, 7] for other invited contributions to the 
CRAS special issue on Dark matter/Dark energy. 

Despite the abundance of new scenarios, one can only really distinguish three 
possibilities: Either the acceleration of the Universe is driven by a Cosmological 
Constant, or dark energy is dynamical, in which case the acceleration of the Universe 
relies on the existence of new degrees of freedom (dof), or finally, the acceleration 
of the Universe relies on dofs already accounted for in Cosmology (e.g. from the 
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backreaction of inhomogeneities). Furthermore, independently of whether or not 
dark energy is dynamical, new degrees of freedom are usually required to tackle the 
Cosmological Constant problem and explain why the vacuum energy is either so much 
smaller than theoretically anticipated or why its gravitational effect is so weak. 

Independently of how these new degrees of freedom tackle the Cosmological Con- 
stant problem, the late-time acceleration of the Universe or possibly the Coincidence 
problem, their mere presence raises a new puzzle of its own: How could they possi- 
bly affect physics on cosmological scales and yet remain hidden on smaller distance 
scales. From energy scales of about lCT 13 eV all the way up to TeV scales, the pres- 
ence of an additional scalar field which would mediate a fifth force is extremely well 
constrained, be it from solar system tests, laboratory tests of General Relativity, or 
cosmological tests such as Big Bang Nucleosynthesis, to name only a few. 

The range of a field is given by its Compton wavelength, thus a scalar field can 
only be relevant on today's cosmological scales if its mass is of the order or smaller 
than the Hubble parameter today. With such a small mass, it is surprising that such 
a degree of freedom would not have already been detected. Over the years, several 
frameworks have been developed to explain how and why light scalar fields could 
be relevant on cosmological scales and yet not manifest themselves on solar system 
scales and every day experiments. The most direct way to hide such a scalar is 
simply to prevent any coupling between that field and the standard model. However 
coupling to gravity usually makes this framework radiatively unstable unless the field 
is protected by a symmetry, such as for instance the discrete shift symmetry present 
in pseudo-Goldstone bosons, which could provide a radiatively stable realization of 
quintessence. 

If, on the other hand, the scalar field couples to matter, the field should be 
hidden via a so-called "screening" mechanism. One of the first realization of such a 
mechanism was proposed by Damour and Polyakov, [8] within the "Least Coupling 
Principle" , providing a relaxation mechanism where the dilaton effectively decouples 
from matter. Besides the Least Coupling Principle, only three screening mechanisms 
are known: The Chameleon [9], the Symmetron [10, 11] and finally the Vainshtein 
mechanism, [12], see Ref. [13] for a review. The underlying physics and screening 
behind these three different mechanisms is very distinct and we shall summarize 
them briefly in this review, however they all have in common the characteristic of 
changing behavior depending on the medium (or on the energy scale of the medium). 
In the Chameleon, for instance, its effective mass depends on the energy density of 
the environment. In the Symmetron and Vainshtein mechanisms, on the other hand, 
it is their effective coupling to matter (or external sources) which is affected by the 
medium, although the explicit realization in the Symmetron and in the Vainshtein 
are very different. 

Whilst any of these previous mechanisms may be realized for scalar fields in their 
own right, the Vainshtein mechanism is also expected to be naturally realized in 
modified theories of gravity where the graviton propagates more degrees of freedom, 
and in particular in theories of massive gravity and their variants. In this review we 
will first summarize in section 2 these three screening mechanisms and then focus 
on the Vainshtein effect. The Vainshtein mechanism is seen explicitly at work in 
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"Galileon" theories and the rest if this review will hence be devoted to Galileons, 
starting by its different realizations in section 3 before moving to its impact on 
Cosmology in section 4. 

2 Screening Mechanisms 

There is little doubts that new physics will appear beyond the Standard Model, some 
of which, such as Dark Matter are expected to have significant effects on the Universe. 
Nevertheless, these new fields are typically expected to be relatively massive and thus 
of short range. These new particles are thus usually expected to bear little effects 
on large cosmological distance scales, of the order of 10 2_3 Mpc or greater. When 
dealing with dark energy, on the other hand, only fields which remain coherent on 
such large scales are expected to have an effect, which is the reason why light scalars 
have been the primary choice. However the strength with which such new degrees 
of freedom can couple to the rest of the standard model is very tightly constrained 
by searches for fifth forces and violations of the weak equivalence principle. These 
tests put constraints on the strength of the interactions mediated by such fields to 
be many orders of magnitude weaker than gravity itself. 

Rather than tuning the coupling between new fields and the Standard Model 
to be small, we take here a different approach and explore the possibility for some 
scalar fields to keep a natural coupling to the Standard Model, whilst still being in 
agreement with observations thanks to a screening mechanism. As we shall see, even 
though only three different classes of mechanisms are known to date, their explicit 
realization can take many facets. However they all share in common the existence of 
non-linear interactions which are essential for the screening to work. 

Starting with a generic scalar field coupled to an external source T (typically the 
trace of the stress-energy tensor of the external sources) , the action for a canonically 
normalized scalar field 4> at the linearized level is* 

£__»<*)»_ i m V + ] £-*r, (1) 

where a is the dimensionless coupling constant. Considering for instance a local- 
ized source at r = 0, T = MS^ 3 \r), then this field typically contributes to the 
gravitational potential with an additional Yukawa-type contribution of the form 
V* = At this level, the only way this field can mediate a force which 

is many orders of magnitude smaller than the gravitational one is either to make it 
very massive (submillimeter tests of Newton's law then impose m >TeV) or to make 
its coupling to matter (unnaturally) small, a <C 1. However we can also imagine a 
situation where the characteristics of the field are medium-dependent, for instance 
the effective mass or the effective coupling (or both) can depend on the energy den- 
sity of the local environment, either locally m = m(p), a = a(p) as will be the case 
in the Chameleon or the Symmetron, or through non-local effects as will be the case 
in the Galileon. 

*We focus on local and Lorentz invariant scenarios. In this review we also focus on single scalar field 
scenarios for simplicity, although most of these frameworks are easily generalizablc to multiple fields. 
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In particular we will be considering a generic scalar field, conformally coupled to 
matter with potential V and wave function Z^, 

C = - X -Z^\4,)d^d v 4, - V(cf>) + A(cf>)T . (2) 

In the three mechanisms presented below, we will see how different screening mech- 
anisms may be realized by appropriate choices of Z, V and A, see Ref. [14]. 

2.1 Chameleon 

The Chameleon is one of the simplest model that exhibits an explicit screening mecha- 
nism, [9]. The idea behind the Chameleon is to start with a light field on cosmological 
scales (in low-energy environments) and increase its (effective) mass in dense regions 
(such as within the Galaxy) so that in dense regions, the force mediated by such a 
field is suppressed by a large Yukawa exponential factor. The explicit realization of 
the Chameleon mechanism relies on self-interactions within the field potential (i.e. 
the field's potential should include a contribution different than the simple mass 
term (ft 2 ) as well as a conformal coupling with external sources, while the field has a 
standard kinetic term Z^ u = g^ u (with g^ v is the space-time metric). 

In the presence of non-relativistic matter, we may split the stress-energy into 
a contribution p from the environment, and a contribution 5T from localized test 
particles, T = —p + ST, such that the potential and the conformal coupling in that 
environment combine to give rise to an effective potential of the form V c r(4>) = 
V(4>) + A((j))p. In the Chameleon, the functions A and V are chosen such that the 
effective potential has a stable minimum at <j) = 4>o(p). Then around that minimum, 
4> = (fro + if, the effective Lagrangian for the scalar field perturbations is of the form 

4t me ieon = ~\(d<P) 2 ~ \ml s { P )4> 2 + a(p)<P5T . (3) 
with the coupling a = A'((f)o(p)) and the effective mass, 

ml s {p) = V"{4>z) + A"{4>,)p. (4) 
such that the Chameleon becomes more massive in denser environments. 

2.2 Symmetron 

An alternative to the Chameleon mechanism has been proposed recently by K. Olive 
and M. Pospelov in [10] as well as by K. Hinterbichler and J. Khoury in [11]. This 
"Symmetron" mechanism shares common features with the Chameleon, such as a 
standard kinetic term Z^ = g^, a conformal coupling to matter and non-linearities 
in its potential, but the realization of the screening relies this time on the weakening 
of the coupling to matter. 



4 



As a specific realization proposed in [11], consider a quadratic coupling to matter 
and a standard 'Higgs-like' potential, 

v{4>) = -Iff + lty* (5) 

m = i + (6) 

such that the effective potential seen by the Symmetron depends on the environment 

VeffW = l(£-S)<? + \w\ (7) 

where here again p is the energy density of the environment. At early times, when the 
Universe is sufficiently hot, or within dense regions of the Universe, where p > (pM) 2 , 
the potential has a unique minimum at <j) = (j)o = 0, which preserves the symmetry 
<f> — > —(ft (or its U(l) analogue if the field is complex). The effective coupling to 
matter seen by fluctuations 4> around that minimum is given by a = A'(0) = 0, 
which explains the weak coupling to external sources. However as the Universe cools 
down, the local energy density drops below the critical density p < p c = (pM) 2 , the 
extremum <p = then becomes a local maximum and the field rolls towards one of the 
two minima located at 4> = <j>o(p) with 0q = (p 2 — p/M 2 )/\ leading to a spontaneous 
breaking of the ^-symmetry. In regions of small density, the effective coupling is 
then non-negligible a = A'((/)q) = (fro/M 2 and the field can have a relevant effect on 
its environment. 



2.3 Vainshtein 

Finally the last class of screening mechanism that has been investigated recently relies 
on non-linearities within the kinetic term Z^(0) rather than within the potential. 
This Vainshtein mechanism was first investigated within the context of Fierz-Pauli 
massive gravity, (see refs. [15]) as a way to evade the vDVZ discontinuity [16], (ex- 
plaining how the additional helicity excitations decouple in the massless limit of the 
graviton), [12]. We will review explicit realizations of the Vainshtein mechanism in 
what follows, but start here by reviewing its main features. Since the Vainshtein 
mechanism does not rely on potential interactions, we may simply set the poten- 
tial V{4>) = and consider a standard and natural conformal coupling to matter, 
A(4>) = 4>/Mp\. The essence of the Vainshtein mechanism then lies in the non- 
standard kinetic term which is symbolically of the form, 

z^{4>) ~ %u + -j^dpM + ^(d^f + ■■■ (8) 

At low-energy, the higher contributions to the kinetic term are negligible, such that 
one recovers a standard kinetic term ~ g^ u and the field couples to external 
sources with a Planck scale coupling. At high energy however, when the background 
configuration for cf) is such that <9 2 0o 3> A 3 , or Z(<j>o) » 1, the higher interactions take 
over and the main contribution to the kinetic term of the fluctuations if around that 
configuration are then of the form Z(4>o)(dip) 2 ~ (<9 2 0o/A 3 ) ra (dip) 2 , with n an integer 
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depending on the details of setup. The properly canonically normalized field is then 
symbolically given by <p = tp/ y/ Z(cf) ), such that the effective coupling to matter 
seen by fluctuations (p around that configuration is given by a = (M P1A /Z(0o)) _1 < 
Mpj 1 . At high energy the field is thus strongly coupled to itself (its interactions are 
important), and becomes weakly coupled to external sources. 

The main issue when dealing with higher order kinetic terms is to ensure that 
no Ostrogradski instabilities arise, [17, 18, 19]. The first explicit realization of the 
Vainshtein mechanism free of such an issue was provided in the context of the DGP 
(Dvali, Gabadadze, Porrati) model [20, 21]. As we shall see below, the decoupling 
limit of DGP model is special in that even though higher derivative interactions can 
be important at a relatively small energy scale A <C Mpi, this does not lead to any 
ghost or standard strong coupling pathologies. These features is shared by a more 
general class of "Galileon" theories, which has been seen to arise in different models 
recently. 

3 Galileons 
3.1 DGP 

The first explicit realization of a soft mass theory of gravity was realized within 
the context of the higher dimensional DGP scenario. Embedded in an infinite five- 
dimensional bulk, the DGP model considers a four-dimensional brane with an in- 
duced Einstein-Hilbert term. By imposing a hierarchy between the four and five- 
dimensional Planck scales (Mpi and M5 respectively), M5 <C Mpi, gravity behaves 
four-dimensional at short distance, up to the cross over scale r c = Mpj/M| = m^ 1 , 
before the five-dimensional effects take over at large distance scales, leading to an 
effectively massive (soft mass) theory of gravity on the four-dimensional brane. 

Since gravity is intrinsically five-dimensional in this scenario, the graviton propa- 
gates five degrees of freedom, one of which behaves a scalar from a four-dimensional 
view point. In the limit where five dimension is "switched off", i.e. M^/Mp\ — > 0, one 
expects to recover four-dimensional standard General Relativity, and this additional 
scalar degree of freedom should hence decouple. This is achieved via the non-linear 
interactions of the extra scalar mode, as originally suggested by Vainshtein. 

The realization of this Vainshtein mechanism is best seen within the so-called 
"decoupling limit" where the helicity-2 and -0 mode of the graviton decouple (which 
we denote respectively as h^ v and (j) in what follows). For this we work at low enough 
energy compared to four-dimensional Planck scale so that the self-interactions of the 
helicity-2 mode may be neglected, but keep a class of non-linearities for the helicity- 
mode alive. This is achieved by taking the limit Mpi,M§ — > 00, while keeping 
the strong coupling scale A = (m 2 Mp\) 1 ^ = M|/Mpi for the helicity-0 mode fixed. 
In this limit, the usual helicity-2 mode of gravity can be treated linearly while the 
interactions for the extra mode (j) remain relevant. The resulting effective action for 
that mode is then, [22] 

= 30D</> - -^(d(p) 2 \3<j)+ -rj—(f>T, (9) 
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leading to the equation of motion, 

+ ^3 ((D0) 2 - (d.d^f) = -jL-T . (10) 

This decoupling limit has several remarkable features: 

1. The effective action is local (contains only a finite number of relevant interac- 
tions). 

2. Even though the interactions involve higher derivative terms (included in □</>), 
they enter in a specific combination such that the equations of motion remain 
second order in derivative and the theory is hence free of the Ostrogradski 
instability. 

3. The equations of motion are invariant under shift and Galileon global trans- 
formations (j) —>■</> + c + v^x^, where c and are constant. This symmetry is 
inherited from five-dimensional Poincare invariance. 

4. As expected, the interactions for (ft remain important up to the much lower 
energy scale A <C M-p\- 

5. This higher interactions rely on "irrelevant operators" from a standard EFT 
point of view, and one cannot in general consider them to be large without 
going beyond the regime of validity of the theory. Here the situation is different 
as these interactions are protected against quantum corrections. This makes it 
possible to set the scale A to be small and to rely on these interactions becoming 
large without needing to include other classes of interactions. This is the essence 
of the "non-renormalization" theorem formulated in Refs. [22] and [23], which 
we describe further in the next subsection. 

To see the Vainshtein mechanism at work, we can focus on spherically symmetric 
static sources, [21]. Taking for instance a localized mass M, for which T = -MS® (r) , 
we can then solve for (f> explicitly, and the force ~ 4>' (r) mediated by this field is 
then given by [23] 

focusing on the '+'- branch, we see that at large distance scales, when r 3> ry, we 
recover the standard Newton square law, <p' Q = M/(Mp\r 2 ), while at short distances 
the interaction term dominates Ocfto ^> A 3 an d the force becomes much weaker com- 
pared to the standard gravitational one, </>o(r) ~ ( 2j ^ p[ ^r) 1 ^ 2 for r <C ry, (this force 
vanishes in the massless limit A — > 0). The cross-over scale between both behaviors 
is given by the "Vainshtein" or "strong-coupling" radius ry = (M /Mp\)Ar l . 

Beyond providing an explicit realization of the Vainshtein mechanism, the decou- 
pling limit of DGP also encodes most of the information on the cosmological behavior 
of this model. In particular the existence of a self-accelerating solution can also be 
seen in this decoupling limit, (corresponding to a the '—'-branch in (11)). 

As mentioned previously, the decoupling limit of DGP has several remarkable 
features which makes it possible to exhibit the Vainshtein mechanism without the 
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presence of any ghost or strong coupling issues (only a finite number of interactions 
remain important and the theory remains under control within the Vainshtein radius) . 
This decoupling limit was recently generalized to the most general set of interactions 
that share the same set of features, and particular corresponds to higher derivative 
interactions that satisfy the Galileon symmetry without leading to any ghost problem 
(no more than two derivatives at the level of the equations of motion), hence called 
the "Galileon" interactions, [24]. Within this "Galileon" framework, one can consider 
a total of four possible interactions, 



C 2 = (d<P) 2 (12) 

£ 3 = W) 2 n</> (13) 

£4 = W) 2 ((D<A) 2 - (<W) 2 ) (14) 

£ 5 = (defy 2 {(U^f -ZU^d^) 2 + 2(8 ^f) (15) 



in addition to a tadpole contribution C± = —<f> which satisfies the accidental sym- 
metry. Since their discovery, the Galileon have been shown to belong to the more 
general class of interactions first introduced by Fairlie et. ai, [25, 26], and can be 
shown to satisfy the recursive relation, 

£ n+1 ~-(50) 2 ^ for n>l. (16) 

It is straightforward to check that there cannot be any further invariant beyond n = 5 
in four dimensions because #^£5 is a total derivative. 

The phenomenology of the Galileon has been explored in great depth in the 
literature, and we will review some of their cosmological implications in what follows, 
but their rich phenomenology relies on the essential fact that they also exhibit a 
Vainshtein mechanism and can hence be screened at short distances or high energy. 
The exact behavior of the field within the strong coupling radius depends on the 
existence of these higher interactions terms, but the suppression of the force is present 
for a wide region of parameter space. 

3.1.1 Non-Renormalization Theorem 

As mentioned previously, the presence of interactions at an energy scale A <C Mp\ 
is essential for the viability of Galileon theories and for the Vainshtein mechanism 
to work. Nevertheless these interactions are typically "irrelevant" from a traditional 
effective field theory sense and the theory is hence not renormalizable. This comes as 
no surprises, since gravity itself is not renormalizable, however we are here required 
to work within a regime where these interactions are important (and even dominate 
within the strong coupling region when d 2 (j) <C A 3 ). With this in mind, one may ex- 
pect that the effective field description goes out of control within the strong coupling 
region and is no longer a valid physical description. 

The existence of the Galileon and shift symmetry goes a long way in preventing 
generic local operators from being generated by loop corrections. However we do still 
expect all local operators which are invariant under shift and Galileon transforma- 
tions to be generated at the quantum level. In particular we expect in principle 
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1. that the Galileon interactions will be themselves renormalized. If correct, this 
would imply that the strong coupling scale A could receive large quantum cor- 
rections <5Aqc ~ M*, (where M* is the cutoff scale, M± » A) which would 
hence destroy the entire Galileon framework, 

2. that quantum corrections will also generate operators of the form (d n (j)) m and 
any combinations of these, with n, m > 2 (which do also respect the Galileon 
symmetry in a trivial way). 

The essence of the non-renormalization theorem formulated in Refs. [22] and [23] 
comes in two parts: 1. The Galileon interactions themselves are not renormalized 
at all (and hence the scale A receives no quantum corrections, or in other words 
the Galileon interactions could be important and radiatively stable) and 2. Higher 
derivative operators of the form (d n 4>) m are generated by quantum corrections, but 
there exists a regime of interest for the theory, for which they are irrelevant, i.e. 
within the strong coupling region, the field itself can take large values, <p ~ A, 
d(f) ~ A 2 , d 2 (f) ~ A 3 , but any further derivative of the field is suppressed, d n 4> <C A n+1 
for any n > 3^. In other words, the Effective Field expansion should be reorganized so 
that operators which do not give equations of motion with more than two derivatives 
(i.e. Galileon interactions) are considered to be large and ought to be treated as the 
relevant operators, while all other interactions (which lead to terms in the equations 
of motion with more than two derivatives) are treated as irrelevant corrections in the 
effective field theory language. 

To understand in more depth the first point related to the non-renormalization of 
the Galileon interactions, let us follow the procedure established in [23]. Having the 
Vainshtein mechanism in mind where the field acquires a large background value, we 
consider perturbations around an arbitrary background configuration (f) = (j) c \ + (p. 
Around that background configuration, the second order Lagrangian for the pertur- 
bations is of the form 

1 

2' 

where the kinetic matrix is symbolically of the form, 

Z~l + c 3 ^- + c 4 A6 +c 5 A9 , (18) 

where 03,4,5 are the coefficients for the 3 rd , 4 th and 5 th Galileon interactions. In order 
to extract the real physical scale at which interactions arise, it is wise to work in 
terms of the canonically normalized field <p = \J~Z ip which acquires the effective mass 

*~5-T' (19) 

The Coleman- Weinberg 1-loop effective action is then of the form 

4 



£^ = --V(^ci)<9 / WV, (17) 



A-ioop = mM<Pc) log ~ log -jfl , (20) 



^Triis is similar to the situation in DBI scalar field models, where the field operator itself and its velocity 
is considered to be large <fi ~ A and dip <~ A 2 , but any higher derivatives are suppressed d n (f> <C A Il+1 for 
n > 2, as shown in section 3.4 
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where M + is the cutoff scale, M* » A. The second part of the non-renormalization 
theorem is now more transparent: even though the field may be large, <p c ~ A, 
d 2 4> c ~ A 3 and hence Z is large, the one-loop effective action only depends on the ratio 
dZ/Z. So once again, quantum corrections are under control as long as dZ <C Z, or 
in other words as long as d 3 (f) c ~ A 4 (and similarly for any higher derivative operator). 
Furthermore, if we pay closer attention to the structure of the counterterms arising 
in the 1-loop effective action, we notice that they all come with at leat one extra 
derivative as compared to the original interactions, dZ ~ d 3 (p. Thus no counterterms 
take the Galileon form, and the Galileon interactions are hence not renormalized (the 
Galileon couplings constants may be tuned to any value and remain radiatively stable, 
such a tuning is hence technically natural) . 

Finally, a last ingredient essential for the quantum stability of Galileon theories is 
the fact that for fluctuations ip on top of the background configuration, interactions 
do not arise at the scale A but rather at the rescaled strong coupling scale A = 
\f~ZA which is much larger than A within the strong coupling region. The higher 
interactions for fluctuations on top of the background configuration are hence much 
smaller than expected and their quantum corrections are therefore suppressed. This is 
the essence of the Vainshtein mechanism itself as presented in section 2.3, explaining 
why the interactions with the fluctuations (p on top of a background configuration 
are suppressed in the vicinity of a dense region of space. 

3.2 Massive Gravity 

As seen previously, the decoupling limit of DGP corresponds to a specific realization 
of the Galileon, which exhibits the Vainshtein mechanism without the presence of 
ghost. It was subsequently established that this does not only apply to soft massive 
gravity such as DGP or its higher dimensional extension, but also to hard mass 
gravity, such as New Massive Gravity in three dimensions, and as we shall see to 
Ghost-free massive gravity in four dimensions. In what follows we show explicitly 
how generalizing the Fierz-Pauli theory of massive gravity leads to a Galileon theory 
in its decoupling limit, and review the procedure underlined in [27, 28]. 
We start by implementing GR with a new potential of the form 

d 4 x v 7 ^ -y 1 (R ~ m 2 U(H^ U , g^j) , (21) 
where H^ v is constructed out of the metric and four Stiickelberg scalar fields 4> a , 

= g„ v - Vabd^du^ , (22) 

such that H^ v represents the fluctuations around flat space-time when taking the 
unitary gauge 4> a = x a , H^ u \ < f >a=x a = g^ u — rj^ u = h^ u /Mp\. As a well-known example, 
the Fierz-Pauli mass term corresponds to the specific choice of potential, such that 
in unitary gauge 

M% iy /=g U F p\ r=x a = --^"v^i hpxKp ~ h^hap) . (23) 
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The Fierz-Pauli mass term is free of any Boulware-Deser ghost at quadratic order, 
however it reappears when working at higher order in perturbations (or equivalently 
when working around a curved background). The existence of this Boulware-Deser 
ghost is most easily seen in terms of the Stiickelberg fields. We start by splitting 
them into the canonically normalized helicity-1 and -0 modes, 

M FX m M P1 m 2 v ; 

In what follows we focus solely on the helicity-2 and -0 modes and set A a = (this 
mode does not get sourced at leading order, and setting it to zero is consistent to all 
orders). The tensor is then expressed as 

12 1 

Hjxv = -y^K v + -j^d^d u (j) - j^r] 01 ? d^d a 4>d v d p <i) , (25) 

where similarly as in DGP, the strong coupling scale A is given by A = (Mpim 2 ) 1 / 3 . 
It is now clear that for a generic potential term Uiji^y) ~ H™ u includes higher 
derivatives in 4> of the form (9 2 ^) n which enter with more than two derivatives 
at the level of the equations of motion. This implies that more than two initial 
conditions ought to be provided. This extra information can be interpreted as a 
physical excitation with the wrong-sign kinetic energy, which in turns implies that the 
associated Hamiltonian of this theory is unbounded from below (at the quantum level 
this signals a ghost while at the classical level this corresponds to the Ostrogradski 
instability), [17, 18, 19]. This is nothing else but the appearance of the Boulware- 
Deser (BD) ghost at the non-linear level. 



3.2.1 Decoupling limit 

To avoid this BD ghost, the potential U ought to be constructed with great care. 
When the BD is present, it can already be observed within the decoupling limit, so 
for simplicity, we will start by constructing a theory which is free of the BD ghost 
within the decoupling limit, and later show how the absence of ghost holds in all 
generality. Similarly as in DGP we take the decoupling limit by sending Mp\ — > oo 
and m — > while keeping the scale A fixed. It will be useful to construct the following 
quantity, [28] 

£"„ = £"„(</, H) = 8% - v 7 ^ - g^H au . (26) 
The tensor K. is built precisely so as to reduce to d^d v 4> in the decoupling limit, 

/C^Imp^oo = ^1^=0 = ^"dM . (27) 

Since self-interactions for <j) arise at the scale A which remains fixed, while the ones 
for the helicity-2 enter at the Planck scale, it is straightforward to see that the leading 
contributions in the decoupling limit can be derived from a simple Taylor expansion, 

m 2 M^^U(H,g) = A 3 ^M P1 W|^ =0 + h^X^ + -±-h 2 ■ ■ , (28) 



11 



with 



= M Pl J-^U(H,g)\ h ^ =0 . (29) 

To ensure that the leading contribution is ghost free, we consider the most general 
following potential 

U(g,H) = _([/C] 2 -[/C 2 ])-a 3 ([/C] 3 -3[/C][/C 2 ] + 2[/C 3 ]) (30) 
-Q 4 ([/C] 4 - 6[/C] 2 [/C 2 ] + 3[/C 2 ] 2 + 8[/C][/C 3 ] - 6[/C 4 ]) , 

where square brackets represent the trace of a tensor, [T] = T^, [T 2 ] = T^T^, etc. . . . 
This potential is specifically built such that its leading contribution in the decoupling 
limit £Y 1/^=0 reduces to a total derivative* . To see the first interactions arising in the 
decoupling limit we then need to work at first order in h^ v . As derived in [27, 28], 
the tensor X^ v associated to the potential (30) is expressed by 

x - l -x<n 1 + 3a * x (2) i 3 ( a 3 + 4a 4 ) (3 ) m) 

^ v ~ A 3 ^ A 6 ^ v A 9 ^ v ' 

with 

X$ = <*V-[<%a- (32) 

x$) = ^i u -[m„-\(m 2 -^ 2 ]Wu (33) 

X$ = ^ - [<^]^ + |(m 2 - [<^ 2 ])^m- - |([^] 3 - 3[<& 2 ][<I>] -h 2[<J> 3 ])^(34) 

with <3? M j, = d^d u 4>. Although this seems to differ from a Galileon model at first sight, 
all the Galileon properties enumerated in section 3.1 are equivalently satisfied here. 
This comes actually as no surprise as the resulting action for the field <p is expressible 
in terms of Galileon interactions after appropriate redefinition of h^ u . Upon the field 
redefinition 

h^u = h» v + 2( H^ ~ ~2A^ 3 d^dz/^ , (35) 

the resulting decoupling Lagrangian for gravity amended with the potential (30) is 
given by 

C = (36) 
+ ^ Q(l + 3« 3 ) 2 - \{a, + 4a 4 )) W ) 2 ([d> 2 ] - [*] 2 ) - S^+^h^ 
' (l + 3« 3 )(a3 + 4a 4 )W) 2 ([f] 3 -3[<l>][<l>] 2 + 2[<I> 3 ]) , 
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^One could also consider an additional linear contribution of the form [K] but this is actually a redundant 
operator, the tadpole can always be reabsorbed by field redefinition of the dynamical metric. 
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where £ is the linearized Einstein tensor, £^uh a/ 3 = Oh^ — d a d^h^ + d^d^h — 

■q^viPih — dadph 01 "). We see in particular that when 03 + 40:4 = we recover precisely 
the Galileon set of interactions (with now a coupling to matter which would be of 
the form 4>T + 1 4 1 ^ 3 <pd u <pT^ u ) , and the Vainshtein mechanism can be seen to work 
for a large range of parameters, even beyond the decoupling limit, [35, 36]. 

In the case where 0-3 + 4«4 7^ 0, the helicity-2 and helicity-0 modes are not gener- 
ically diagonalizable in a covariant way. As we shall see below, this has interesting 
consequences for cosmology and represents a new coupling which is not encoded by 
the Galileon family of interactions. 

3.2.2 Beyond the decoupling limit 

The previous subsection focused on how to construct a theory of massive gravity 
which is ghost-free in the decoupling limit, which may only be used to describe 
physical processes that happen at most at the energy scale A. Since this scale is very 
small, (typically A ~ lCT 40 Mpi), one really ought to understand the consistency of 
the theory beyond this decoupling limit. Fortunately we have derived the potential 
for the graviton in a completely covariant way, and so the action (21) for massive 
gravity, implemented with the potential (30), where H^ u is given in terms of the four 
Stiickelberg fields as in (22) remains a fully valid description beyond the decoupling 
limit (the only difference when working beyond the decoupling limit is that higher 
order interactions with h^ v may no longer be omitted). 

To explore the phenomenology of this theory beyond the decoupling limit (i.e. at 
energy scales above A), one can of course keep the same language as what was used 
within the decoupling limit and split the degrees of freedom into two modes present 
in the metric gw = rj^ v + h^ u , and the other modes included in the Stiickelberg 
fields 4> a = x a + V a . Even though valid, this language turns out not to give the best 
description of the system. The reason for this is that the quantity 4> that appears in 
the Stiickelberg field as V a = d a 4> does no longer fully describe the helicity-0 mode 
of the graviton when interactions with gravity are included, [29, 30]. 

To be more precise, when ignoring gravity (i.e. when working around flat space- 
time, which is effectively what happens within the decoupling limit since all higher 
interactions in h are suppressed in the limit Mpi — > 00), space-time enjoys a global 
Lorentz invariance x a — > A^x b . Moreover, the theory is also invariance under an 
independent internal transformation of the four Stiickelberg fields, which is also a 
global Lorentz transformation (f) a — > A^(p b . Thus within the decoupling limit, one 
may work in the representation of the single group A = A, and the Stiickelberg 
fields hence behave as vectors under this global Lorentz symmetry. It does therefore 
makes sense to perform a scalar-vector-tensor decomposition of the fields in this 
representation, and express <fi a = x a + A a + d a <p in that limit, where <p then captures 
the physics of the helicity-0 mode. However beyond the decoupling limit space-time 
loses the global Lorentz symmetry and 4> no longer captures the relevant degree of 
freedom. 

When maintaining the split (f> a = x a + A a + d a <fi beyond the decoupling limit, then 
(f) will maintain its Galileon symmetry in a manifest way (simply by construction), 
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but this will not be the case for the helicity-0 mode itself. Furthermore beyond the 
decoupling limit, <f> may appear with higher derivatives in its equations of motion, 
however once again, since <ft bears no explicit physical meaning beyond the decoupling 
limit (no physical processes involve only scatterings of <j) with itself), this does not 
imply the presence of ghost. Indeed, one can show explicitly that unitarity is still 
preserved, [29, 30]. However to see this more explicitly, it is wiser to work in the 
Hamiltonian language, and to work in unitary gauge 4> a = x a so that all the degrees 
of freedom are then carried by the metric. This analysis was successfully performed 
to fourth in perturbations in the Hamiltonian in Ref. [28], and then generalized to 
all orders in [32] . The existence of a secondary constraint was then demonstrated in 
[33] . Recent work in Ref. [34] also clarifies how the absence of ghost in the decoupling 
limit generalizes beyond that limit. 



3.3 New Massive Gravity 

Whilst a massless spin-2 field propagates no dynamical degree of freedom in three 
dimensions, the same does not hold for a massive spin-2 field which can excite two 
polarizations, one of which behaves as a helicity-0 mode. This realization has made 
the study of three-dimensional massive gravity especially appealing. One explicit 
realization of massive gravity in three dimensions was proposed a few years ago by in 
[37] and relies on the introduction of higher curvature terms. Such terms usually lead 
to new dofs with the opposite sign kinetic term and would be unacceptable in most 
cases, however since massless spin-2 field do not propagate any physical dof, one can 
"survive" in three dimensions with a wrong-sign massless spin-2 field without further 
ado. Whilst the situation is very specific to three dimensions and cannot be extended 
to higher dimensions in an straightforward manner, the resulting decoupling limit 
theory carries very similar features to what is obtained in four-dimensional theories 
of massive gravity. To see this explicitly, we start with the original formulation of 
NMG, 

•Snmg = J d 3 xM P1 ^ (-R + -L (j^ _ l R *)j . (37) 

To see the Galileon structure arising, let us focus on the conformal mode, g^ v = 

^1 + y /w~ ) ^ anc ^ cons ider the decoupling limit Mp\ — > oo, whilst keeping the 

scale A 3 d = (m 2 y/Wp[) 2 ^ fixed. In that limit, the resulting action for NMG is then 
precisely of a Galileon form, [38] 

= 2 » " "4/2 W) 2 ^ . (38) 

Thus in all known models of massive gravity be it DGP, New Massive Gravity or 
the Ghost-free theory of massive gravity in four-dimensions, the absence of the BD 
ghost, the decoupling limit resembles that of a Galileon theory. Whilst these models 
represent the natural embedding of Galileon scalar fields beyond the decoupling limit, 
we shall see in what follows another class of braneworld configuration from where 
Galileons arise when considering a non-relativistic limit. 
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3.4 DBI - Galileon 



Some of the features of the Galileon are very similar in spirit to the DBI world action 
for a brane embedded within a wrapped extra dimension. In particular, DBI inflation 
relies on a braneworld action of the form 



where dimension-8 and other higher irrelevant operators such (dcfi) 4 are important, 
d<j) ;$ /(</ ) ) 2 without yet going beyond the regime of validity of the theory. The reason 
for that is here again due to the existence of a non-renormalization theorem, whereby 
a class of irrelevant operators may be important at energy scales much lower than 
the Planck scale d(f> ~ f 2 <C Mpj but involve quantum corrections which themselves 
are negligible as long as d 2 <p <C f 3 . These features is very similar to the Galileon, 
and as we shall see below, both theories can actually be "reunited" within a common 
braneworld framework as shown in Ref. [39]. 

To see this explicitly, let us consider a four-dimensional brane, embedded within 
a five-dimensional maximally symmetric extra dimension (we focus on flat extra 
dimension within this review, although this framework is easily generalizable to more 
general geometries). The Lovelock invariants have already been established almost 
half a century ago [40] as providing all the possible geometrical quantities which are 
free of any ghostly additional degree of freedom. In four dimensions, there are two 
of them, namely the cosmological constant and the Einstein-Hilbert scalar curvature 
term, leading to the usual Einstein equation. We can therefore start by considering 
these two contributions on the four dimensional brane. In five dimensions, in addition 
to these two geometrical contributions, the Gauss-Bonnet term represent a non-trivial 
higher-derivative term which combines the curvature terms in such a way so as to 
remain second order at the level of the equations of motion. If present within the 
fifth dimension, the Ricci scalar curvature and the Gauss-Bonnet term induce specific 
boundary terms on the brane similarly as the Gibbon-Hawking boundary term. In 
addition to these four geometrical contributions, we can also include a tadpole term 
arising from the integration over the bulk loops, similarly as in DBI. Considering all 
these geometrical quantities, the most general four-dimensional action arising from 
all the possible five and four-dimensional Lovelock invariant is thus of the form 



where £2 is the brane tension or four-dimensional cosmological constant which we 
choose to write in terms of a scale m of dimension mass and a contribution from 
the tadpole, and £3,4,5 are respectively the contributions for the extrinsic curvature, 
the four-dimensional Einstein-Hilbert term and the boundary term associated with 




(39) 



(40) 



S = M|[ / & A xyf^g (£ 2 + c 3 £ 3 + c 4 £ 4 + c 5 £ 5 ) , 




(41) 
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Gauss-Bonnet curvature in the bulk, 

C 2 = m 2 (^=-^J (42) 

£ 3 = mK (43) 
£4 = R (44) 
£5 = m- l K GB . (45) 

We consider a flat bulk with a four-dimensional brane localised at y = (f> = <f)/M-p\m. 
The induced metric on the brane is then given g^ u = 7]^ u + d^(pd u (j), leading to the 
extrinsic curvature on the brane, 

K, u = --Md=, (46) 

V 1 + (W 

whilst the boundary term associated with the Gauss-Bonnet curvature in the bulk is 
[41]: 

/Cgb = -\k%, + KK% - l -K* - 2G^ V , (47) 

where G^ u is the Einstein tensor associated with the induced metric g^ v . The exact 
expression for these different invariants was provided in Ref. [39]. Since the brane 
is embedded within a flat fifth dimension, the four-dimensional action should be 
invariant under translation and boosts along the fifth dimension, which applied on 
the brane position, corresponds to the infinitesimal global symmetry, 

<f> -» (ft + c + v^x" + ^ 2M 2 . (48) 

We now consider the non-relativistic limit, (dcf>) — > 0, which can be achieved, by 
sending once again Mpi — > 00 while keeping the scale A = (Mpim 2 ) 1 / 3 fixed. The 
resulting action on the brane is then nothing else than the four Galileon interactions 
proposed in (12-15), with the associated non-relativistic global Galileon symmetry: 

4> -)• 4> + c + . (49) 

A analogue construction of the Galileon interactions from five-dimensional Lovelock 
invariants was also later developed using a compactified extra dimension a la Kaluza- 
Klein, [42]. 

In all of these models, be it braneworld models or hard or soft massive gravity, 
the scalar new degree of freedom which encodes either the helicity-0 mode or the 
brane-bending mode behaves as a Galileon scalar field within a specific "decoupling" 
limit. When expressed in terms of a Galileon scalar field, the Vainshtein mechanism 
is most easily tractable, and all these models enjoy a healthy strong coupling effect, 
at least for some choice of their parameters. This strong coupling is essential in 
understanding how such a light scalar mode could be present and remain as yet 
unobserved. In what follows we consider different cosmological applications of the 
Galileon and show how it can play an important role for the late-time acceleration 
of the Universe as well as for inflation. 
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4 Galileon Cosmology 



4.1 Self-acceleration in DGP and Galileons 

One of the breakthrough of the DGP model, [20], was the realization that the Universe 
could be accelerating without the presence of any cosmological constant or other 
source of "dark-energy" but sourced by the graviton very own dofs and in particular 
by its helicity-0 mode, [43]. Within this "self-accelerating" solution of DGP, the 
expansion of the Universe is naturally accelerating without the need of any external 
source or cosmological constant. It was quickly realized though that this branch of 
solution contained a negative norm state which was hence propagating a ghost. To 
see this explicitly, let us consider the decoupling limit of DGP which takes the form 
(including the helicity-2 mode), 

4 d GP = -\h a ^^h, u + ± ± (d<pfn^ + l-V^ + ^ T > (50) 

where h^ u and 4> represents the helicity-2 and -0 modes respectively, and the linearized 
induced metric on the brane in DGP is given by g^ v = rj^ + (h^ + 2^^) . The 
'—'-sign corresponds to the normal (stable) branch of DGP while the '+'-branch is 
the self-accelerating one. This leads to the standard Einstein's equations for the 
helicity-2, 

M Pl SG^ = £$h af) = j^-T^ , (51) 

such that in the absence of any source the helicity-2 mode is not excited (h^ u = 
up to gravitational waves). For the helicity-0 mode on the other hand, there is a 
non-trivial configuration even in the absence of external sources, 

3D0T ^ ((D0) 2 - (W) 2 ) = -J^ T - ( 52 ) 

The homogeneous configuration d 2 <p = is always a solution when T = 0, but 
there exists also another branch which admits 4>q = ±^A 3 x fl x tJ ' as a solution. For 
this branch of solution, the linearized induced metric on the brane is of the form 
9uu = (1 i A^x 2 )^^, which is nothing else but (anti) de Sitter linearized around 
Minkowski. In particular when working with the ' — '-sign, one can only achieve a 
Minkowski or an AdS solution, while the '+'-sign allows for a de Sitter solution. 
This goes to show how the induced metric on a DGP brane can be de Sitter even 
in the absence of any cosmological constant or dark-energy. However one can easily 
notice that fluctuations on the top of this self-accelerating solution bear no kinetic 
term, already signaling the presence of strong coupling issues, on top of the self- 
accelerating branch. This issue is resolved when including any external source, be it 
radiation, dark matter or other typically present in the Universe, but it can then be 
shown that the helicity-1 mode of the graviton then bears a negative kinetic term 
hence signaling the existence of a ghost, [44]. 

The presence of a ghost within the self-accelerating branch of DGP can be avoided 
when working in a more general Galileon framework, developed in Ref. [24] (or other 
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extensions of DGP, [45]). The equation of motion for the helicity-0 mode (52) then 
gets replaced by 

Cl U<t> + g (m 2 - [$ 2 ]) + ^ ([Of - 3[<!> 2 ][<I>] + 2[$ 3 ]) 

+^(m i -mV]+^ 2 } 2 +8[m 3 }-^ 4 }) = -^r, (53) 

so there are non-trivial branches of solution in the absence of matter T = 0, for 
different choices of parameters 02,3,4. In particular there exists a region in parame- 
ter space for which fluctuations around the self-accelerating solution are stable and 
exhibit the Vainshtein mechanism, [24]. 

In what follows, we focus instead on a slightly different realization of self-accelerating 
solutions which relies on the new specific interactions that arise within massive grav- 
ity. 



4.2 The accelerating Universe in Massive Gravity 

To study the self-accelerating solutions in massive gravity, we review the approach 
provided in [46]. In this case it is more convenient to work in Jordan frame where 
the linearized metric is given by g^ v = rj^ + h^/M-pi and h^ v is sourced both by 
matter and the helicity-0 mode, 

M Pl 5G, u = £$h af) = j±- (T, v + TW) , (54) 

where T^ v is the additional contribution arising from the helicity-0 mode, 

tM = A 3 A>. (55) 
The equation of motion for <p is given by: 

flWfcV^ ("^V + ^^V"Vr + 3(Qt8 2 t, 4fl4) ^^) = 0> (56) 

where e^ VOL p is the Levi-Cevita tensor. This has two important consequences, first 
of all the existence of an additional T^J implies the existence of self-accelerating 
solutions in the absence of matter, T^ u = 0, as along as T^' ~ g^ u . However, 
another peculiarity of this framework is the existence of degravitating or screening 
solutions, where the helicity-0 mode can absorb the contribution from a cosmological 
constant, T^ v = —\cc9txv and T^ u = AccflWj i n such a way that the geometry 
remains flat (g^ u = 77^) even in the presence of a potentially large cosmological 
constant. Such a solution would not be possible within a pure Galileon framework 
and the difference originates from the additional coupling to matter d fl (f)d u (j)T^ L ' that 
arises when switching to the Einstein frame. As we shall see below such terms can 
also have important observational signatures. 
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4.2.1 Self-accelerating Branch 

We start by giving a brief overview of the self-accelerating solution, before moving 
onto the degravitating branch. We focus on linearized solutions around flat space- 
time and look for configurations of the form h^ u = —^H 2 x 2 r] lll/ , where H is the 
equivalent of the Hubble parameter and <p = ^A 3 qx 2 . This branch corresponds to 
the case where the constant q is related to the Hubble parameter by (56) 

H 2 + 2a 2 q + 3a 3 q 2 ^ = , (57) 

with a 2 = — 5(1 + 3^3) and 03 = §(03 + 4a4), and the Einstein equation in the 
absence of external sources is given by (54) 

M Pl H 2 = 2qA 3 (- i + a 2 q + a 3 q 2 ^j . (58) 

Self-accelerating solutions for which H > thus exist for a large family of parame- 
ters 02,3. To investigate this branch further, we analyze small perturbations for the 
helicity-2 and -0 modes, h^ v = — \H 2 x 2 "q lJLl , + x^v and <p = \A?qx 2 + (p. The resulting 
Lagrangian for the perturbations (up to a total derivative) reads as follows [46] 

C = -^X^^Xa/3 + 6(a 2 + 3a 3 g)^^^-3a3^^(^) 2 n^ 

+ ^^x-4l ) M + ^x-4i ) M + ^. (59) 

As mentioned in [46], a key feature of this self-accelerating branch is the fact that 
there is no quadratic mixing between x an d <p, which implies that x identifies the 
correctly diagonalized helicity-2 mode at quadratic order and only the helicity-2 mode 
couples to external sources T^ v at that level. Thus on top of the self-accelerating 
solution, the helicity-0 mode does not couple to external sources at quadratic order 
and will only do so at higher order (through a coupling of the form d fl ipd u (pT^' / ). This 
results remains unaffected by the presence of cosmological matter at the background 
level. Therefore, for arbitrary external sources, it is always consistent within the 
decoupling limit to keep the fluctuation of the helicity-0 unexcited, ip = 0. This 
mechanism is a unique feature of the self-accelerating branch in Ghost-free massive 
gravity. Beyond the decoupling limit we expect the coupling between the helicity-0 
and -2 mode to re-emerge (and hence a coupling of the helicity-0 mode to matter), 
however such a coupling ought to be suppressed by additional powers of m, and 
should therefore only be relevant a very large distance scales. 

From the perturbed Lagrangian (59), we see directly that the situation is quite 
different from DGP. In particular even though the helicity-2 and -0 modes do not 
couple at the linearized level, the helicity-0 mode keeps a non-trivial kinetic term 
proportional to a 2 + ?>a 3 q. Fluctuations around the self-accelerating branch are hence 
healthy (no ghost nor tachyon) if a 2 + 3a 3 q > 0. This condition, combined with the 
background equations of motion (57) and (58) and the requirement that H 2 > are 
satisfied simultaneously if the parameters 02,3 satisfy [46] 

a 2 < and - ^ < a 3 < . (60) 
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In particular we see that the presence of the third tensor X^J (which has direct equiv- 
alent in pure Galileon theories) is essential for the stability of this self-accelerating 
branch. The resulting Hubble parameter is then of the order of the graviton mass, 

H 2 = m 2 [2a 2q 2 + 2a 3 q 3 - q }>0, with q = + ^ + . (61) 

003 3y2a3 

Further investigations of this self-accelerating branch have established that the helicity- 
1 mode is likely to become unhealthy around this configuration unless it remains 
unexcited, [47]. This would have serious implications for the self-accelerating branch 
which should be understood further. We now turn to the other branch of solution 
which allows the degravitation of the cosmological constant. 



4.2.2 Degravitating branch 

An orthogonal approach to self-acceleration, is to address the question of whether 
a large cosmological constant (vacuum energy) could be present, whilst only given 
rise to a very weak late-time acceleration of the Universe. This is the idea behind 
degravitation, [48]. A first step in the study of degravitation is to understand whether 
there can exist a static solution (Minkowski space-time) even in the presence of a 
cosmological constant, which could be a late-time attractor for the Universe. One 
way to achieve this within massive gravity is to "absorb" the effect of the cosmological 
constant through the helicity-0 mode. We review here again the framework derived in 
[46], and consider a cosmological constant, T^ v = —Xccd^u absorbed by the helicity-0 

(it) 

mode T^u = Xqc9^u in such a way that the resulting curvature vanishes, = 0. 
Using the same notation as for the self-accelerating branch, we hence set H = for 
which (57) is then trivially satisfied and the helicity-0 mode satisfies 

1 , 2 , 3 ^cc rao s 

~2 q + a2q +asq =-6A^- (62) 

As long as the parameter 03 is present (0:3 + Aa^ / 0), Eq. (62) has always at least 
one real root. There is therefore a flat solution for an arbitrarily large cosmological 
constant. A further analysis of the scalar, vector and tensor fluctuations shows that 
this branch of solution if it exists is always stable. 

However a drawback of this degravitating solution is that the resulting scale for 
helicity-0 interactions are no longer governed by the parameter A, but rather by the 
scale determined by the cosmological constant to be absorbed, A ~ (Acc/^pi) 1 ^ 3 - 
To see this, let us pursue the analysis of the fluctuations around the degravitating 
branch and keep the higher order interactions. The resulting Lagrangian involving 
the helicity-0 mode is then 



with 



£ {2) = - V [<p] + |x£> M + m) , (63) 

a 2 a 2 + 3a 3 q M P1 



A 3 A 3 {-l+Aa 2 q + 6a 3 q 2 y A CC 

a 3 2a 3 / Mpi 



A6 A 6 (-l + 4a 2 g + 6a 3 g 2 ) 3 V^cc 



(64) 
(65) 
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assuming 02,3 ~ 0(1). Fifth force constraints usually requires interactions to kick in 
at least at a scale A 3 /Mpi < (1CP 33 eV) 2 (this bound might be relaxed slightly in 
the presence of X^), which implies that the the maximum allowed value of vacuum 
energy that can be screened without being in conflict with observations is fairly low, 
of the order of (1CT 3 eV) 4 , (unless the coefficients 02,3 are tuned to different natural 
values, which would still remain technically natural). Although not phenomenolog- 
ically viable, the existence of this degravitating solution opens the door for a new 
set of possibilities to tackle the cosmological constant problem. This degravitating 
solution can also see an analogue within the full theory beyond the decoupling limit, 
[49]. See Refs. [50] for further work in massive gravity with an open FRW cosmology, 
Ref. [51] for self-accelerating solutions in the full massive gravity theory, Refs. [52] 
for cosmology in bi-gravity theories and Refs. [53] for cosmology in massive gravity 
beyond the decoupling limit. 

4.3 Galileon Inflation 

As seen in the previous section, the Galileon can play an essential role for cosmology 
at a time where the Hubble parameter is of the order of H ~ (A 3 /Mpi) 1 / 2 , such that 
for very low strong coupling scale A ~ lCT 13 eV the Galileon affects the late evolution 
of the Universe, whereas if the Galileon interactions arise at a much larger strong 
scale, they could have an important impact during inflation (in that case the Galileon 
should then decay during reheating explaining why it is not observed today, since 
a larger strong coupling scale would not allow for a sufficiently strong Vainshtein 
mechanism at late time). This philosophy has received a great attention in the past 
few years, and one cannot within this short review do justice to all the advances 
made in that direction. We therefore only focus on the model of Galileon Inflation 
presented in [54], and refer to [55] for further work on Galileon Inflation, G-Inflation 
and the generation of non-gaussianities in these classes of models and to Refs. [56, 57] 
for work on Galilean Genesis. 

The idea behind Galileon Inflation is to rely on the non-renormalization theorem 
of the Galileon class of interactions to build a radiatively stable model of inflation. 
We hence consider the Galileon to play the role of the usual inflaton scalar field, 
with all the Galileon interactions (or their covariant counterpart around an FRW 
background, [58, 59, 39]) arising at the scale H <C A <C Mp\ (where H is the Hubble 
parameter during inflation) in addition to the tadpole contribution X(f>, which satisfies 
the same properties as the other Galileon interactions. All these operators satisfy 
the shift symmetry which is essential for the production of a scale invariant power 
spectrum, but as such would not allow for a graceful exit of inflation. However in 
addition to these operators one may also consider a pure mass term m 2 (j) 2 which 
breaks the shift symmetry but does affect the non-renormalization theorem in any 
way, [54] . The stability of the effective field theory under radiative corrections is thus 
preserved in the presence of this mass term, which now allows for a graceful exit of 
inflation. 

These considerations are technically correct only in flat space-time while the cou- 
pling with gravity breaks the Galilean symmetry, which in turn leads to additional 



21 



corrections. However these corrections are suppressed by at least two powers of 
{H/ A) <C 1 and can thus be neglected during inflation. 

In the case where the Galileon interactions are negligible Z = He/)/ A 3 <C 1, 
the perturbations are weakly coupled and one recovers a standard slow-roll inflation 
scenario. However Galileon Inflation also allows for another regime Z > 1, where the 
interactions can be important without yet going beyond the regime of validity of the 
Effective Field Theory. This is the regime we are interested in, and the cosmological 
background evolution is then determined by the largest Galileon operator (assuming 
all the Galileon interactions arise at the same coupling scale A). 

Focusing the study on the effect of fluctuations around this cosmological back- 
ground we follow small fluctuations £, constructed from the background solution by 
mapping t i— > t + £(x, t) on a hypersurface of constant time and work in a gauge 
where such hypersurfaces are spatially flat. Up to cubic order in these fluctuations, 
the resulting Lagrangian then includes terms of the form, 



£ € 5 a{tf 



a(t) 2 v s ' J /s a(t) 2SV s ' a{tf 



, (66) 



where a{t) is the scale factor and 31,2,3 as well as a(t) and c s (t) depend on the 
background field configuration and can vary over a relatively large range of scales 
depending on the exact value of the dimensionless scale Z. The operators with 
coupling g\ and 32 are dimension-6 operators, while the last one with coupling g% is 
a dimension-7 operator. In a standard field theory, as well as for DBI inflation, one 
cannot rely on the dimension-7 operator to be important without going beyond the 
regime of validity of the Effective Field Theory. In the case of Galileon Inflation, 
the situation is however different, and the operator-7 can well be of the same order 
as the others if the dimension-6 operators are suppressed by additional powers of 
H/A. In the case where the dimension-7 operator is non-negligible, the structure of 
the non-gaussianities is much richer than what is present in other standard models 
of slow-roll inflation or even DBI inflation. In particular while the two dimension- 
6 operators £ 3 and £(<9£) 2 typically lead to a contribution to the non-gaussianity 
parameter /jy^ ~ gi/a and similarly the operator /jy^ ~ g2/ac 2 s respectively, the 
dimension-7 operator is suppressed by even more powers of the sound speed, /jy^ ~ 
gs/ac^. DBI inflation relies on the second contribution which is already suppressed 
by two powers of the sound speed. In the case of Galileon Inflation, the amount of 
non-gaussianity could be of the same order or even higher, but carried by a different 
set of operators and hence leading to different shapes of triangle. 

Whilst the exact predictions of this model are quite sensitive to the precise value of 
the different parameters, it opens up the possibility of a new class of non-gaussianities. 
This possibility was then further investigated by [60] if the leading interactions from 
operators of the form (d<j)) 2 {d 2 <j)) n vanish, the next to leading order operators of the 
form (d 2 <fi) n become the dominant ones. Dimension-9 operators can then lead to 
important contributions to the non-gaussianities and arise with their own specific 
shape. More precisely, two of these operators can lead to equilateral triangles whilst 
another operator leads to a flattened isocele triangle. 
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5 Outlook 



Infrared modified theories of gravity typically come with new degrees of freedom 
which manifest themselves at low-energy. At high energy, on the other hand we 
expect to recover a standard theory of gravity and these new dofs are generically 
screened via a Vainshtein mechanism. This mechanism can be seen explicitly at work 
in braneworld models or models of massive gravity, such as DGP, NMG in three 
dimensions or Ghost-free massive gravity in four-dimensions. In all these models, 
the helicity-0 mode resembles in some limit a Galileon scalar field and can have an 
interesting phenomenology for the late-time acceleration of the Universe, be it self- 
acceleration or degravitation, if the strong coupling scale of the Galileon interactions 
is very low. For parametrically larger values of the strong coupling scale, the Galileon 
could on the other hand affect the very early Universe and in particular be a candidate 
for the inflaton, leading to non-trivial non-gaussianities in the power spectrum of 
curvature fluctuations. 

The screening of the helicity-0 mode at high energy implies that the departure 
from General Relativity on solar system scales are very much suppressed. Neverthe- 
less, this screening can yet have distinct signatures in cosmology and in particular 
for structure formation, [61]. Another key signature of such models is the deviation 
from the standard Friedmann equation, that could distinguish these models from 
other scenarios of modified gravity or with additional dynamical degrees of freedom, 
see for instance Refs. [62, 63]. 

Another effect which can lead to specific signatures in models of massive gravity, 
is the emergence of a new type of coupling to matter of the form d /1 (j)d u <fiT llu . Whilst 
this coupling is not relevant for spherically symmetric static sources as it vanishes in 
that case, it can be significant for relativistic sources. It has already been established 
by M. Wyman in Refs. [64] and [65], that its effect in weak lensing could provide a 
smoking gun signature for this class of models. 

Finally one peculiarity of Galileon models is the generic appearance of superlu- 
minal propagation, at least when the Vainshtein mechanism is at work, which may 
lead to the creation of Close-Timelike Curves (CTCs), [66, 67]. For instance if one 
of the dimensions were compact, CTCs may easily be created without the need of 
any additional exotic type of matter, however whenever one reaches a regime where 
CTCs may form, the Galileon inevitably becomes infinitely strongly coupled which 
leads to the breakdown of the effective field theory used to described the formation 
of CTCs in the first place. In addition, the formation of CTCs requires a background 
solution which is unstable with an arbitrarily fast decay time, similarly to a ghost- 
like instability. This suggests that an analogue of Hawking's chronology protection 
conjecture is also in effect for Galileons, [68]. 

References 

[1] S. M. Carroll, Living Rev. Rel. 4, 1 (2001) [arXiv:astro-ph/0004075]. 

[2] E. J. Copeland, M. Sami and S. Tsujikawa, Int. J. Mod. Phys. D 15, 1753 (2006) 
[arXiv:hep-th/0603057]; 



23 



[3] R. Durrer and R. Maartens, Gen. Rel. Grav. 40, 301 (2008) [arXiv:0711.0077 
[astro-ph]]. 

[4] P. Astier and R. Pain, "Observational Evidence of the Accelerated Expansion 
of the Universe" , Invited contribution for the CRAS special issue on Dark mat- 
ter/Dark energy. 

[5] M. Kunz, "The Phenomenological Approach to Modeling Dark Energy", Invited 
contribution for the CRAS special issue on Dark matter/Dark energy. 

[6] J. Martin, "Everything You always Wanted to Know about the Cosmological 
Constant (but Were Afraid to Ask)", Invited contribution for the CRAS special 
issue on Dark matter /Dark energy. 

[7] C. Clarkson, 'Establishing Homogeneity of the Universe in the Shadow of Dark 
Energy", Invited contribution for the CRAS special issue on Dark matter/Dark 
energy. 

[8] T. Damour and A. M. Polyakov, Nucl. Phys. B 423, 532 (1994) [hep- 
th/9401069]. 

[9] J. Khoury, A. Weltman, Phys. Rev. Lett. 93, 171104 (2004). [astro-ph/0309300]. 

[10] K. A. Olive and M. Pospelov, Phys. Rev. D 77, 043524 (2008) [arXiv:0709.3825 
[hep-ph]]. 

[11] K. Hinterbichler, J. Khoury, Phys. Rev. Lett. 104, 231301 (2010). 
[arXiv: 1001.4525 [hep-th]]. 

[12] A. I. Vainshtein, Phys. Lett. B 39, 393 (1972). 

[13] J. Khoury, arXiv:1011.5909 [astro-ph. CO]. 

[14] A. J. Tolley, "Dark Energy", Cosmo '09. 

[15] M. Fierz and W. Pauli, Proc. Roy. Soc. Lond. A 173, 211 (1939); W. Pauli and 
M. Fierz, Helv. Phys. Acta 12, 297 (1939). 

[16] H. van Dam and M. J. G. Veltman, Nucl. Phys. B 22, 397 (1970); V. I. Zakharov, 
[Pisma Zh. Eksp. Teor. Fiz. 12, 447 (1970)]. 

[17] M. Ostrogradski: Mem. Ac. St. Petersbourg VI 4, 385 (1850) 

[18] E.T. Whittaker, fourth ed, Cambridge University Press, Cambridge, 1937, pp. 
265267. 

[19] R. P. Woodard, Lect. Notes Phys. 720, 403 (2007). 

[20] G. R. Dvali, G. Gabadadze and M. Porrati, Phys. Lett. B 485, 208 (2000) 
[arXiv:hep-th/0005016]; G. R. Dvali and G. Gabadadze, Phys. Rev. D 63, 
065007 (2001) [arXiv:hep-th/0008054]. 

[21] C. Deffayet, G. R. Dvali, G. Gabadadze and A. I. Vainshtein, Phys. Rev. D 65, 
044026 (2002) [arXiv:hep-th/0106001]. 

[22] M. A. Luty, M. Porrati, R. Rattazzi, JHEP 0309, 029 (2003). [arXiv:hep- 
th/0303116 [hep-th]]. 

[23] A. Nicolis, R. Rattazzi, JHEP 0406, 059 (2004). [hep-th/0404159]. 



24 



[24] A. Nicolis, R. Rattazzi, E. Trincherini, Phys. Rev. D79, 064036 (2009). 
[arXiv:0811.2197 [hep-th]]. 

[25] D. B. Fairlie, J. Govaerts, A. Morozov, Nucl. Phys. B373, 214-232 (1992). [hep- 
th/9110022]; D. B. Fairlie, J. Govaerts, J. Math. Phys. 33, 3543-3566 (1992). 
[hep-th/9204074]; D. B. Fairlie, J. Nonlin. Math. Phys. 9, 256-261 (2002). [math- 
ph/0112033]. 

[26] D. Fairlie, [arXiv: 1102. 1594 [hep-th]]. 

[27] C. de Rham, G. Gabadadze, Phys. Rev. D82, 044020 (2010). [arXiv:1007.0443 
[hep-th]]. 

[28] C. de Rham, G. Gabadadze, A. J. Tolley, Phys. Rev. Lett. 106, 231101 (2011). 
[arXiv: 101 1.1232 [hep-th]]. 

[29] C. de Rham, G. Gabadadze and A. Tolley, arXiv: 1107.3820 [hep-th]. 

[30] C. de Rham, G. Gabadadze and A. J. Tolley, JHEP 1111, 093 (2011) 
[arXiv:1108.4521 [hep-th]]. 

[31] S. F. Hassan, A. Schmidt-May and M. von Strauss, arXiv:1203.5283 [hep-th]. 

[32] S. F. Hassan and R. A. Rosen, arXiv:1106.3344 [hep-th]. 

[33] S. F. Hassan and R. A. Rosen, arXiv:1111.2070 [hep-th]. 

[34] M. Mirbabayi, arXiv:1112.1435 [hep-th]. 

[35] K. Koyama, G. Niz and G. Tasinato, Phys. Rev. D 84, 064033 (2011) 
[arXiv:1104.2143 [hep-th]]. 

[36] G. Chkareuli and D. Pirtskhalava, arXiv: 1105. 1783 [hep-th]. 

[37] E. A. Bergshoeff, O. Hohm, P. K. Townsend, Phys. Rev. Lett. 102, 201301 
(2009). [arXiv:0901.1766 [hep-th]]. 

[38] C. de Rham, G. Gabadadze, D. Pirtskhalava, A. J. Tolley and I. Yavin, JHEP 
1106, 028 (2011) [arXiv: 1103. 1351 [hep-th]]. 

[39] C. de Rham and A. J. Tolley, JCAP 1005, 015 (2010) [arXiv:1003.5917 [hep-th]]. 

[40] D. Lovelock, J. Math. Phys. 12, 498 (1971). 

[41] S. C. Davis, Phys. Rev. D 67, 024030 (2003) [arXiv:hep-th/0208205]. 

[42] K. Van Acoleyen and J. Van Doorsselaere, Phys. Rev. D 83, 084025 (2011) 
[arXiv: 1102.0487 [gr-qc]]. 

[43] C. Deffayet, Phys. Lett. B 502, 199 (2001) [arXiv:hep-th/0010186]; C. Deffayet, 
G. R. Dvali and G. Gabadadze, Phys. Rev. D 65, 044023 (2002) [arXiv:astro- 
ph/0105068]. 

[44] K. Koyama, Phys. Rev. D 72, 123511 (2005) [hep-th/0503191]; D. Gorbunov, 
K. Koyama and S. Sibiryakov, Phys. Rev. D 73, 044016 (2006) [hep-th/05 12097]; 
C. Charmousis, R. Gregory, N. Kaloper and A. Padilla, JHEP 0610 (2006) 066 
[hep-th/0604086]. 

[45] C. de Rham and A. J. Tolley, JCAP 0607 (2006) 004 [hep-th/0605122]. 



25 



[46] C. de Rham, G. Gabadadze, L. Heisenberg and D. Pirtskhalava, Phys. Rev. D 
83 (2011) 103516 [arXiv:1010.1780 [hep-th]]. 

[47] K. Koyama, G. Niz and G. Tasinato, JHEP 1112, 065 (2011) [arXiv:1110.2618 
[hep-th]]. 

[48] G. Dvali, G. Gabadadze and M. Shifman, Published in *Minneapolis 2002, 
Continuous advances in QCD* 566-581; arXiv:hep-th/0208096; G. Dvali, 
G. Gabadadze and M. Shifman, Phys. Rev. D 67, 044020 (2003) [arXiv:hep- 
th/0202174]; N. Arkani-Hamed, S. Dimopoulos, G. Dvali and G. Gabadadze, 
arXiv:hep-th/0209227; G. Dvali, S. Hofmann and J. Khoury, Phys. Rev. D 76, 
084006 (2007) [arXiv:hep-th/0703027]. 

[49] G. D'Amico, C. de Rham, S. Dubovsky, G. Gabadadze, D. Pirtskhalava and 
A. J. Tolley, Phys. Rev. D 84, 124046 (2011) [arXiv: 1108.5231 [hep-th]]. 

[50] A. H. Chamseddine and M. S. Volkov; Phys. Lett. B 704, 652 (2011) 
[arXiv: 1107.5504 [hep-th]]. A. E. Gumrukcuoglu, C. Lin and S. Mukohyama, 
JCAP 1111, 030 (2011) [arXiv:1109.3845 [hep-th]]; A. E. Gumrukcuoglu, C. Lin 
and S. Mukohyama, arXiv: 1111.4107 [hep-th]. 

[51] K. Koyama, G. Niz and G. Tasinato, Phys. Rev. Lett. 107, 131101 (2011) 
[arXiv:1103.4708 [hep-th]]. 

[52] S. F. Hassan and R. A. Rosen, rXiv:1109.3515 [hep-th]; M. S. Volkov, JHEP 
1201, 035 (2012) [arXiv:1110.6153 [hep-th]]. M. von Strauss, A. Schmidt-May, 
J. Enander, E. Mortsell and S. F. Hassan, arXiv: 111 1.1655 [gr-qc]; D. Comelli, 
M. Crisostomi, F. Nesti and L. Pilo, arXiv: 111 1.1983 [hep-th]. 

[53] C. de Rham and L. Heisenberg, Phys. Rev. D 84, 043503 (2011) [arXiv:1106.3312 
[hep-th]]. 

[54] C. Burrage, C. de Rham, D. Seery and A. J. Tolley, JCAP 1101, 014 (2011) 
[arXiv: 1009.2497 [hep-th]]. 

[55] T. Kobayashi, M. Yamaguchi and J. 'i. Yokoyama, Phys. Rev. Lett. 105, 
231302 (2010) [arXiv: 1008.0603 [hep-th]]; S. Mizuno and K. Koyama, Phys. 
Rev. D 82, 103518 (2010) [arXiv: 1009.0677 [hep-th]]; A. De Felice, R. Kase and 
S. Tsujikawa, Phys. Rev. D 83, 043515 (2011) [arXiv: 101 1.6132 [astro-ph.CO]]; 
A. De Felice and S. Tsujikawa, JCAP 1104, 029 (2011) [arXiv:1103.1172 
[astro-ph.CO]]; T. Kobayashi, M. Yamaguchi and J. 'i. Yokoyama, Phys. 
Rev. D 83, 103524 (2011) [arXiv:1103.1740 [hep-th]]; Z. -G. Liu, J. Zhang 
and Y. -S. Piao, Phys. Rev. D 84, 063508 (2011) [arXiv: 1105.5713 [astro- 
ph.CO]]. T. Kobayashi, M. Yamaguchi and J. 'i. Yokoyama, Prog. Theor. 
Phys. 126, 511 (2011) [arXiv:1105.5723 [hep-th]]; X. Gao, JCAP 1110, 021 
(2011) [arXiv:1106.0292 [astro-ph.CO]]; X. Gao and D. A. Steer, JCAP 1112, 
019 (2011) [arXiv:1107.2642 [astro-ph.CO]]; A. De Felice and S. Tsujikawa, 
Phys. Rev. D 84, 083504 (2011) [arXiv:1107.3917 [gr-qc]]; S. Renaux-Petel, 
Class. Quant. Grav. 28, 182001 (2011) [Erratum-ibid. 28, 249601 (2011)] 
[arXiv:1105.6366 [astro-ph.CO]]; S. Renaux-Petel, S. Mizuno and K. Koyama, 
JCAP 1111, 042 (2011) [arXiv:1108.0305 [astro-ph.CO]]; X. Gao, T. Kobayashi, 



26 



M. Yamaguchi and J. 'i. Yokoyama, Phys. Rev. Lett. 107, 211301 (2011) 
[arXiv:1108.3513 [astro-ph.CO]]; 

[56] P. Creminelli, A. Nicolis and E. Trincherini, JCAP 1011, 021 (2010) 
[arXiv: 1007.0027 [hep-th]]. 

[57] L. Levasseur Perreault, R. Brandenberger and A. -C. Davis, Phys. Rev. D 84, 
103512 (2011) [arXiv:1105.5649 [astro-ph.CO]]; T. Qiu, J. Evslin, Y. -F. Cai, 
M. Li and X. Zhang, JCAP 1110, 036 (2011) [arXiv:1108.0593 [hep-th]]. 

[58] C. Deffayet, G. Esposito-Farese and A. Vikman, Phys. Rev. D 79, 084003 (2009) 
[arXiv:0901.1314 [hep-th]]. 

[59] C. Deffayet, S. Deser and G. Esposito-Farese, Phys. Rev. D 80, 064015 (2009) 
[arXiv:0906.1967 [gr-qc]]. 

[60] P. Creminelli, G. D'Amico, M. Musso, J. Norena and E. Trincherini, JCAP 
1102, 006 (2011) [arXiv:1011.3004 [hep-th]]. 

[61] N. Chow, J. Khoury, Phys. Rev. D80, 024037 (2009). [arXiv:0905.1325 [hep-th]]; 

[62] N. Afshordi, G. Geshnizjani and J. Khoury, JCAP 0908, 030 (2009) 
[arXiv: 08 1 2 . 2244 [astro-ph] ] . 

[63] B. Jain and J. Khoury, Annals Phys. 325, 1479 (2010) [arXiv: 1004.3294 [astro- 
ph.CO]]. 

[64] M. Wyman, Phys. Rev. Lett. 106, 201102 (2011). [arXiv: 1101. 1295 [astro- 
ph.CO]]. 

[65] S. Sjors and E. Mortsell, arXiv:1111.5961 [gr-qc]. 

[66] A. Adams, N. Arkani-Hamed, S. Dubovsky, A. Nicolis and R. Rattazzi, JHEP 
0610, 014 (2006) [hep-th/0602178]. 

[67] J. Evslin and T. Qiu, JHEP 1111, 032 (2011) [arXiv:1106.0570 [hep-th]]. 

[68] C. Burrage, C. de Rham, L. Heisenberg and A. J. Tolley, arXiv: 11 11.5549 [hep- 
th]. 



27 



